Self-similar relaxation dynamics of a fluid wedge in a Hele-Shaw cell 
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Let the interface between two immiscible fluids in a Hele-Shaw cell have, at t = 0, a wedge 
shape. As a wedge is scale-free, the fluid relaxation dynamics are self-similar. We find the dynamic 
exponent of this self-similar flow and show that the interface shape is given by the solution of an 
unusual inverse problem of potential theory. We solve this problem analytically for an almost flat 
wedge, and numerically otherwise. The wedge solution is useful for analysis of pinch-off singularities. 
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Introduction. Interface dynamics between two immis- 
cible fluids in a Hele-Shaw cell have attracted a great 
interest in the last two decades. Most of the efforts have 
dealt with forced flows, when a more viscous fluid is dis- 
placed by a less viscous fluid. In the forced case the 
viscous fingering instability Q, develops and brings 
about intricate issues of pattern selection in a channel 
geometry 0, [3, development of fractal structure in a 
radial geometry ,6j , etc. The role of small surface tension 
in the theory of a forced Hele-Shaw flow is to introduce 
regularization on small scales. This Letter deals with an 
unforced Hele-Shaw (UHS) flow 0, H [H EH , where 
surface tension at the fluid interface is the only driv- 
ing factor. The pertinent free boundary problem here is 
non-integrable and, because of its non-locality, hard for 
analysis. To our knowledge, the only known analytical 
solutions to this class of problems are (i) a linear analy- 
sis of the dynamics of a slightly deformed flat or circular 
interface 0,0 and (ii) a recent asymptotic scaling anal- 
ysis of the dynamics of a l ong stripe of an inviscid fluid 
trapped in a viscous fluid To get more insight into 
the physics of UHS flows, we address here the case when 
one of the fluids at t = has the form of a wedge. In 
this case the flow is self-similar. Building on this simpli- 
fication, we recast the problem into an unusual inverse 
problem of potential theory. We solve this problem an- 
alytically for an almost flat wedge and numerically for 
several other wedge angles. Finally, we use a wedge so- 
lution for analysis of pinch-off events of the UHS flow, 
which has attracted much interest in theory and experi- 
ment 00. 

Governing equations and self- similarity. Let one of the 
fluids have a negligible viscosity, so that the pressure in- 
side this fluid is constant and can be taken zero. The ve- 
locity of the viscous fluid is v (r, t) = — (& 2 /12/i) Vp (r, t), 
where p is the pressure, /i is the dynamic viscosity, and 
b is the plate spacing 0,0 0- Therefore, the interface 
speed is 

v n = -(b 2 /l2fi)d n p, (1) 

where index n denotes the components of the vectors 
normal to the interface outwards, and d n p is evaluated 
at the respective points of the interface 7. As V • v = 



y 




FIG. 1: The setting for fluid wedge relaxation. 

in the (incompressible) viscous fluid, the pressure there 
is a harmonic function: 

V 2 p = . (2) 

The Gibbs-Thomson relation at the interface yields 

p| 7 = (7r/4)<7/C, (3) 

where a is surface tension, and K. is the local curvature of 
the interface, positive when the inviscid region is convex 
outwards. As the flow is undriven we demand 

d n p = at r^oo. (4) 

We assume that the interface has the form of a graph 
y = y{x, t) and rewrite Eq. Q as an evolution equation: 

d t y(x,t) = -(b 2 /12^)d nP ^l + (d x y) 2 

= (6 2 /12/i) [d x y(x, t) d xP - d y p] , (5) 

where the derivatives of p are evaluated at the interface. 
At t = the inviscid fluid has the form of a wedge of 
angle a, so that y — — \x\ cot (a/2), see Fig. ^ As this 
initial condition and Eqs. 101- © do not introduce any 
length scale, the solution must be self-similar [12j. Let 
L(t) be the retreat distance of the wedge tip. Then the 
interface position and the pressure in the viscous fluid 



2 



can be written as 
y(x,t) 
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where primes stand for A"-derivatives. Now, using Eqs. 
© and Q in Eq. JSJ), we arrive at the following equation: 

&(X)d x P-dYP = \[$(X)-X&(X)], (9) 

where the derivatives of P are evaluated at the rescaled 
interface $(X), = [(7rAcrfe 2 i)/(16/i)] 1/3 , and A is an 
unknown dimensionless parameter. The boundary con- 
ditions are $(0) = -1, $'(0) = and <f>(X -> ±oo) = 
— \X\ cot (a/2). Note that we have already found the dy- 
namic scaling exponent 1/3: the same as observed in the 
relaxation of fractal viscous fingering patterns 0, Ho| . 
The shape function &(X) (and the parameter A) for a 
given wedge angle a is determined by the solution of 
the following (quite unusual) inverse problem of poten- 
tial theory. A harmonic function P(X, Y) must obey both 
a Dirichlct boundary condition [Eq. JSJ], and a Neu- 
mann boundary condition [Eq. iJSJ], while the function 
$(X) must be determined from the demand that these 
two conditions be consistent. We solved this problem an- 
alytically for an almost flat wedge, and numerically oth- 
erwise. Before reporting the analytic solution, we present 
a large-X asymptote of &(X), valid for any wedge an- 
gle. It corresponds to the leading term of the multipole 
expansion of P(X, Y) at large distances. Introduce, for 
a moment, polar coordinates r, <f> with the origin at the 
point X = Y = and measure the polar angle <f> from the 
ray Y — —Xcot (a/2) counterclockwise. At large \X\ the 
curve Y = $(X) is almost fiat, so P[X, $(X)] -> there 
by virtue of Eq. (JSJ) . Therefore, the leading term of the 
multipole expansion is P(r 3> 1,0) = const r~ u sin(i/0), 
where v = (2 — a/n)^ 1 13]. Now we employ Eq. © and 
obtain, at \X\ ^> 1, 



$(A) = -\X\ cot (a/2) + C \X\- 



(10) 



with an unknown constant C that depends only on a. 

Almost flat wedge. Let us assume that 7r — a <C vr, 
and introduce the small parameter e = cot(a/2) <C 1. 
We rescale the variables: X = £/e, F = rj/e, P(X, F) = 
e 2 U(^,n) and A = Ae 3 . In the rescaled varaibles the 
interface equation is -q = £i/>(£), where ipiO = &(£./ £ )- 



Keeping only leading terms, we can rewrite the boundary 
conditions (JSJ and © for the harmonic function U(£,rj) 
in the following form: 



respectively. We fix the coordinates by choosing y(x — where 
0,t) = —L(t), that is $(0) = —1. In the rescaled coordi- 
nates X — x/L(t) and Y — y/L(t) the Laplace's equation 
© keeps its form, while Eq. Q becomes 



U[t,e1>(0] = ^'(0, (11) 

WfceiKO] = akv'(0-^(0] . (i 2 ) 



1&(0) = -1, V'(0) = 0, 4>(£ -» ±oo) = -|£| + o(l) . (13) 



The rescaled problem does not include e, except in the 
second argument of the functions on the left hand side 
of Eqs. l|ll[l and (|12fl . In view of the condition — > 
±cx)) = — 1£| one cannot put the second argument to zero 
at sufficiently large |£|. As will be shown below, these 
values of £ are exponentially large in e -1 , while at shorter 
distances one can safely put the second argument to zero. 

The problem obtained in this way is soluble exactly. 
Assume ^>(£) is known. Then one can easily find the 
harmonic function in the upper half-plane r/ > 0, that 
satisfies the Dirichlet condition it(£,0) = 011 the 

£-axis: 



1 f°° r]ip"(s)ds 

* J -co (i - s) 2 + v 2 



(14) 



Now we should impose the Neumann condition 112|) 
(where we put e = 0). To avoid calculation of hyper- 
singular integrals, we find the harmonic conjugate 
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(15) 



and, by virtue of the Cauchy-Riemann conditions, replace 
d n U(£, 0) by — d^V(^, 0). This yields a non-standard 
integro-differential equation 

W®- m --l'f (16, 

where ^ denotes the principal value of the integral. For- 
tunately, upon differentiation with respect to £ Eq. 1|16|) 
becomes an equation for if>" (£) which is soluble by Fourier 
transform. The result is 
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(the constant of integration is determined from the con- 
dition J_ t/j"(£)d£ = —2). Integrating twice in £ and 
using the first two conditions in Eq. (|13|) yields 



m = -i - - 



• 2 fcf 

sin 



(18) 



To determine A, we expand this expression at |£| 3> 1 
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where T(w) is the gamma- function |l4|. To elimi- 
nate the offset 0(1) we put A = (8/3) 7r" 3 [r(2/3)] 3 = 
0.213545 . . .. Though the integrals in Eqs. |JT7J| and 
(|T%)) can be expressed via the generalized hypergeomet- 
ric function ri F n (a; b; z), it is more convenient to keep the 
integral form [15j. To complete the solution, we find the 
rescaled pressure: 

U{£,v) = / e-^r-^ cos k£, dk. (20) 

n J-oo 

Now we find the distance |£| = 1(e) ^> 1 at which the 
solution (|18|) becomes inaccurate, and improve the large- 
ly asymptote. Let us compare Eq. |JTSJ, which becomes 

^) = -iei + 4[r( ^ 2 /3)] 3 r 2 , i« i^i «/(£), (21) 

with the large- 1^| multipole asymptote (|10|l : 

Tp(0 = -\t\ + C(e)\{;/e\-^ , |e|»l- (22) 

where, for small e, — (3ir — a)/(2ir — a) ~ — 2 + 2e/tt. 
We see that the last term in Eq. (|21|l lacks the small 
correction 2e/n in the exponent of £. We can match the 
two asymptotes (12111 and ill'l'l) in their common region of 
validity 1< |£| < /(e). We define 1(e) as the value of |£| 
for which the correction to the exponent yields a factor e: 
1(e) = e"/' 2 ^ [notice that, at |£| ~ 1(e), the deviation of 
ip(£) from its flat asymptote — 1£| is already exponentially 
small: ~ e _7r / e ]. The matching yields C(e), and we arrive 
at the improved small-e large- 1£| asymptote: 

^-^ + www e ' ¥ '"' r ' t¥ ' (23) 

So far we have dealt with inviscid fluid wedges: a < 
180°. Our results, however, can be immediately extended 
to viscous fluid wedges: a > 180°. 

Numerical algorithm and parameters. For a general 
wedge angle the shape function of the self-similar inter- 
face can be found numerically. Instead of dealing with 
the similarity formulation of the problem (JSJ-JHl, we com- 
puted the time-dependent relaxation of wedges of differ- 
ent angles, as described by (rescaled) Eqs. 0-@ [lif . 
Our numerical algorithm |l7| employs a variant of the 
boundary integral method for an exterior Dirichlct prob- 
lem for a singly connected domain, and explicit tracking 
of the contour nodes. The harmonic potential is repre- 
sented as a superposition of potentials created by a dipole 
distribution with an unknown density D on the contour. 
D is computed from a linear integral equation [is| . Com- 
puting another integral of this dipole density yields the 
harmonic conjugate, whose derivative along the contour 
is equal to the normal velocity of the interface. 

We chose the singly connected domain to be (i) a rhom- 
bus with angles 120° and 60°, (ii) a square, and (iii) a 




FIG. 2: The retreat distance L(t) and its power-law fit for 
a = 120° (triangles), 90° (squares) and 60° (circles). 



straight cross with aspect ratio 10 3 |16|. In this man- 
ner we could exploit the 4-fold symmetry of the domains 
and measure the retreat distance of the respective ver- 
texes, L(t), and the rescaled interface shapes &(X) for 
four wedge angles: 120°, 90°, 60° and 270°, the latter one 
corresponding to a 90° wedge of the viscous fluid. The 
ultimate shapes of the rhombus- and square-shaped do- 
mains are perfect circles. Therefore, to observe the self- 
similar asymptotics we did the measurements at times 
much shorter than the characteristic time of relaxation 
toward a circle, and at distances much smaller than the 
domain size (so that the effect of the other vertexes could 
be neglected). For the rhombus and square an equidis- 
tant grid with 901 nodes per side was employed. For 
the quarter of the cross we used 2801 nodes. The time 
step was taken to be 10 -3 times the maximum of the 
ratio of the interface curvature radius and the interface 
speed at the same node. The domain area conservation 
was used for accuracy control. For the measurements 
reported here the area was conserved with an accuracy 
better than 10 _3 %. 

Numerical results. We first report the results for the 
three viscous fluid wedges. Figure shows the retreat 
distance L(t) for the angles 120°, 90°, and 60°. Power 
law fits yield L(t) = 0.48 i 33 , 0.84 t 033 , and 1.33 * 33 , 
respectively, so the dynamic exponent 1/3 is clearly ob- 
served. In the rescaled units, used in the simulations 
[lr| . the analytical prediction for an almost flat wedge 
is L(t) = at 1 / 3 , where a = [3A) 1/3 e ~ 0.862 e. For 
a = 120° and 90° this yields a ~ 0.498 and a = 0.862, 
respectively, in very good agreement with the measured 
values 0.48 and 0.84. Even for a = 60° the analytical pre- 
diction, a — 1.493, is only 12% higher than the measured 
value 1.33. 

The rescaled shapes of the three evolving wedges are 
depicted in Fig. That the curves, measured at three 
different times, collapse into a single curve proves self- 
similarity. The prediction of our almost-flat-wedge the- 
ory, shown on the same three graphs, works very well for 
a = 120° and 90°, and fairly well even for 60°. 

We also measured, for each of the three values of the 
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FIG. 3: (color online). Right panel: the shape function $ 
for wedges of inviscid fluid: a = 120° (a), 90° (b) and 60° 
(c). Data for three different times (2.3 x 10" 7 , 1.3 x 10" 5 
and 1.1 x 10" 4 for a and c, and 8.2 x 10" 7 , 8.9 x 10" 6 and 
3.5 x 10~ 5 for b) collapse into a single curve. The red solid 
line is the prediction of the almost-flat- wedge theory, the blue 
dashed line is the asymptote Y = — \X\ cot a/2. Left panel: 
a blowup of a part of figure c. 



shape functions of the two viscous fluid wedges. The case 
of (3 = 90° is shown in Fig.01 where the retreat distance, 
the shape function and the pinch-off time are taken from 
the previously described simulation of the cross-shaped 
domain with aspect ratio 10 3 . 

Summary. We have studied analytically and numeri- 
cally the surface tension driven flow of a fluid wedge in 
a Hele-Shaw cell. We have shown that the fluid inter- 
face evolves self-similarly, found the asymptotic interface 
shape at large distances, and recast the problem into an 
unusual inverse problem of potential theory. We solved 
this inverse problem analytically in the limit of nearly 
flat wedge, and performed numerical simulation which 
support and extend the analytic calculations. Like in the 
case of self-similar solutions, obtained for wedge-like ini- 
tial conditions in other surface tension driven flows |l9j |. 
this solution provides a sharp characterization of the UHS 
flow. It also sheds a new light on the pinch-off singulari- 
ties of this flow. 




FIG. 4: (color online). Pinch-off of a straight branch of thick- 
ness A of inviscid fluid. Each of the two viscous fluid wedges 
corresponds to j3 = 90° (i.e. a — 270°). The blue dashed 
lines: the interface shape at t = 0; the red solid lines: the 
interface shape at the pinch-off time t = 50.65 (in the units of 
[T^I L The inset: the measured viscous fluid retreat distance 
versus time and its power-law fit 0.74 x t ' 33 . 



wedge angle, the tail of the shape function (the difference 
between $(X) and Y = —\X\ cota/2). The results are 
in excellent agreement with the theoretical prediction, 
given by the last term in Eq. i|l(J|) . 

Pinch-offs. The self-similar wedge solutions are useful 
for analysis of pinch-offs in UHS flows 0, Q . Let the in- 
viscid fluid domain represent, at t = 0, an infinitely long 
straight branch, coming at an angle f3 from an infinitely 
long straight "trunk" . The simple physics in the inviscid 
fluid branch precludes interaction between the two vis- 
cous fluid wedges of angles [3 and tt — (3, which evolve in a 
self-similar manner, causing the inviscid branch to thin, 
and ultimately to pinch-off. The i 1 / 3 law intrinsic in the 
self-similar solution implies that the pinch-off time is pro- 
portional to the branch thickness cubed. The interface 
shape at all times prior to the pinch-off can be obtained, 
with a proper rescaling, from the respective self-similar 
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